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ABSTRACT 

Recent analyses of the fluctuations of the soft Diffuse X-ray Background (DXB) have 
^ I provided indirect detection of a component consistent with the elusive Warm Hot Inter- 

galactic Medium (WHIM). In this work wc use theoretical predictions obtained from 
hydrodynamical simulations to investigate the angular correlation properties of the 
WHIM in emission and assess the possibility of indirect detection with next-generation 
X-ray missions. Our results indicate that the angular correlation signal of the WHIM 
is generally weak but dominates the angular correlation function of the DXB outside 
ly^ , virializcd regions. Its indirect detection is possible but requires rather long exposure 

ly^ ■ times [0.1-1] Ms, large (~ 1° x 1°) fields of view and accurate subtraction of isotropic 

fore/background contributions, mostly contributed by Galactic emission. The angular 
correlation function of the WHIM, which turns out to be positive for 9 < 5' provides 
limited information on its spatial distribution. A satisfactory characterization of the 
WHIM in 3D can be obtained through spatially resolved spectroscopy. 1 Ms long ex- 
posures with next generation detectors will allow to detect ~ 400 O vii-|-0 vni X-ray 
emission systems that could be used to trace the spatial distribution of the WHIM. We 
k> , predict that these observations will allow to estimate the WHIM correlation function 

j_j ■ with high statistical significance out to ^ 10 Mpc h~^ and characterize its dynamical 

' state through the analysis of redshift-space distortions. The detectable WHIM, which 

is typically associated with the outskirts of virialized regions rather than the filaments 
has a non-zero correlation function with slope 7 = — 1.7±0.1 and correlation length 
To = 4.0 ± 0.1 Mpc h~^ in the range r = [4.5, 12] Mpc h~^. Redshift space distances 
can be measured to assess the dynamical properties of the gas, that we predict to be 
typically infalling onto large virialized structures. 

Key words: diffuse radiation - large-scale structure of universe (WHIM) - X-rays: 
diffuse background 
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1 INTRODUCTION 

The observed abundance of baryons in the local Uni- 
verse does not agree with the high redshift observa- 
tions (|Fukugita. Hogan. fc PeeblesI 119981 ). The Lyman al- 
pha forest at 2 = 2, the results from the WMAP 
experiment, and the predictions of the standard nucle- 
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osynthesis model are all consistent with a baryonic den- 
sity D,B ~ 0.04 jRauch et al.l [l99l: [Weinberg et al.l Il997 



Buries fc Tvtleij Il998l: iKirkman et al.l l2003l: iBennett et al l 



20031 : ISpergel et al.ll2007l : iKomatsu et al.ll2009l ). At variance 
with these results, the baryon mass in stars, galaxies, and 
clusters in the local Universe seems to be 2 — 4 times lower. 

Large-scale cosmo l ogical hydrodyn a mic simulations 
(ICen fc Ostrikerl[l999al: ICroft et al.ll200ll: IVoshikawa et"al 
I2OO3I : IChen et albooj iBorgani et al.ll2004llcen fc Ostrikei 
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20061: IShen. Wadslev fc StinsonI I2OI0I : ISmith et alj I2OI0I : 
Tornatore et all I2OI0I )" allow to follow the evolution of 



baryons from high redshift to present days. A common fea- 
ture of these simulations is that 40 — 50 % of the baryons 
in the Universe at 2 ~ are in the form of highly ionized 
gas with temperature in the range of 10^ — 10^ K. This 
gas is characterized by a spatial distribution that closely 
resembles that of the underlying cosmic web of the Dark 
Matter, and is commonly indicated as the Warm-Hot Inter- 
galactic Medium (WHIM). According t o thermal emission 
models for gas in coUisional equilibrium (|Ravmond fc SmithI 
[l923) and assuming an average metallicity of 0.1 Z©, this 
highly ionized gas emits mainly in the soft X-ray ([0.1- 
1] keV) and in the Far and Extreme Ultra Violet (FUV, 
[0.01-0.1] keV) energy bands. However, the expected WHIM 
contribution to the total surface brightness (SB) in these 
bands is expected to be weak. The comparison between the 
predicted emission of the WHIM and the measured Dif- 
fuse X-Ray Background (DXB), for example, shows that 
the WH IM contribution to the total SB is in the order 
of 10% (iPhillips, Cen, fc Ostrikcr"200l': ' Ursino fc Galeazzil 
I2OO6I : |u rsino. Roncarelli. fc Galcazzi 2 010| ). Highly ionized 
metals in the warm-hot gas are responsible for line emis- 



pected 


to be rather weak dKravtsov, Klvoin. fc Hoffman 


2002 




Chen ct al. 20031: Klvoin et al. 2003: Viel et al. 


2003 


Yoshikawa et al.1 |2003|. |2004|: ICen fc Fand l2006l: 


Branchini et al.ll2009l). 



Observational data confirm the difficulty to detect 
WHIM lines, both in absorption and emission. A clear 
detection of the WHIM comes from the absorption lines 
of ions like O vi, N v, C iii, C iv, Si iii. Si iv, and 
Fe III in the FUV spectra of bright, distant sources. 
Roughly a hundred FUV absorbers have been detec ted in 
the spectra of ~30 AGN (Daufprth fc S hull 200^, I2OO8I : 
iThom fc Cheiil l2008l : iTripp et al.1 120081 ). Their abundance 
turne d out to be in ag reement with theoretical predic- 
tions (|Cen fc Fand I2OO6I ) . However, these absorbers trace 
only the warm gas (T< 10*' K) whereas the bulk of the 
WHIM is expected to be at higher temperatures. In the 
soft X-ray band t he search for the WHIM has l ed to 



FauE. Canizares. fc Yao 20071: Buote et al. 20091: Fang et al.1 


2OIOI: WiUiams et al. 


2010|) whose statistical significance is 


still under debate ( 


Kaastra et al.l I2OO6I: iReismussen et al.l 



In emission, the only claim for a line detection so far 
is likely to be associated with a gas fil ament connecting th e 
two galaxy clusters A222 and A223 (|Werner et al.ll2008l ). 
However, the overdensity (5 ~ 150) and the temperature 
(T~ 1 X 10^ K) of this gas are only mar ginally consistent 
with t he WHIM defini t ion. I n addition, Zappacosta et ahl 
l|2002h : [Mannucci et all (|2007^ observed a soft X-ray excess 
associated to galaxy concentration and possibly consistent 
with WHIM emission. 

Finally, an indirect evidence for the WHIM comes 
from the analysis of the Diffuse DXB in the soft band 
[0.4-1.0] keV, which is known to be mainly contributed 
by unresolved AGN, starburst galax ie s and local dif- 
fuse e mission (ISnowden et aL Il997 ; Mushotzkv et al] 



2OOOI: iMcCammon et al.l |2002|: iGaleazzi et al.1 
Henley fc Shehon I2OO8I : iGupta fc Galeazzil 
Henlev fc ShehonI I2OI0I) . These 
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sources account for all 



but 10% to 20% of the DXB (jHickox fc MarkevitchllioOTh 
matching the theore tical expectations fo r the WHIM 
emission in this band (jRoncarelli et al.]|2006l ). 

The WHIM is also expected to provide some contri- 
bution, perhaps with characteristic features, to the angu- 
lar correlation properties of the DXB fluctuations. Such 
contribution, if detected, would provide a further, indirect 
evidence for the WHIM. Indeed, a recent angular correla- 
tion analysis of the unresolved DXB measured by XMM 
(jGaleazzi. Gupta, fc Ursindl2009l ) has revealed the presence 
of a component that has been identified with the WHIM. Fu- 
ture X-ray missions are expected to detect a large number 
of characteristic emission lines in the X-ray band through 
which one will be a ble to unambiguously detect the WHIM 
(jCen fc Fandl2006l ) and probe its spatial distribution much 
more densely than it could be possibly done in absorption, 
due t o the paucity of the background sources ijViel et al.l 
l2003l , l2005l ). Therefore one expects that only emission stud- 
ies in the X-ray band will be able not only to detect the 
WHIM and characterize its thermal state, but also to pro- 
vide a tomography of a large fraction of the missing baryons 
in the local universe. In particular, as we will show in this 
work, given the characteristic of next-generation X-ray de- 
tectors, the estimate of the spatial two-point correlation 
function of the WHIM line-emitting regions will allow to 
characterize the spatial correlation properties of the WHIM 
and its dynamical state. 

In this work we investigate the possibility of charac- 
terizing the angular and spatial distribution of the WHIM 
using next-generation X-ray satellites. More specifically, we 
aim at two difi'erent but related goals. First we investigate 
the possibility of finding the signature of the WHIM in the 
angular correlation function of the unresolved DXB. Sec- 
ond, we assess the possibility of tracing the spatial distri- 
bution and the dynamical state of the WHIM by detecting 
its characteristic emission lines and measuring the spatial 
two-point correlation function of the line-emitting regions. 
The two ingredients required for this goal are: a theoreti- 
cal WHIM model to rely upon and a reference experimental 
se tup. For the WHIM m odel we assume t he sa me one as 
of lUrsino et"al] (I2OIOI') an d lBranchini et al.l (|2009l '). based on 
the Borgani et al. (200J) hydrodynamical simulation. As a 
reference experimental setup we consider the recentl y pro- 
posed JfeniaQ X-ray satellite (jHartmann et al]|2009l ). The 
reason for considering this particular mission is that it is de- 
signed to carry two instruments well suited for our purposes: 
the High Angular Resolution Imager (HARI) and the Cryo- 
genic Imaging Spectrometer (CRIS). The HARI is a CCD or 
CMOS based detector with a field of view (f.o.v) shaped as a 
circle of diameter 1.4°, angular resolution of 15" for on-axis 
objects, and energy resolution of 50 — 150 eV FWHM (in 
the range 0.3 — 5.0 keV). The CRIS is a spatially resolved 
spectrometer based on Transition Edge Sensor technology. 
The imaging area is 0.9° x 0.9°, angular resolution is ~ 2.5', 
and the energy resolution is in the range [1-3] eV. The X- 
ray telescope has an effective area of 1000 cm^ at 1 keV. 
The HARI instrument is expected to provide the possibility 
to study the angular correlation properties of the DXB and 
isolate the WHIM contribution thank to its large f.o.v. and 
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angular resolution, smaller t han the typical angular s ize of 
the WHIM emitting element (|Ursino fc GaleazzJlioO^ ') . The 
CRTS instrument will allow to detect the WHIM emission 
lines and identify the spatial position of the emitting gas 
thanks to the good angular and energy resolution. 

The plan of the paper is as follows. In Section [2] we 
introduce the WHIM model and the procedure we adopted 
to generate the simulated X-ray maps and spectra that will 
be considered for the correlation analysis. In Section [3] we 
describe the estimator used to compute the angular two- 
point correlation function and estimate the contribution of 
the WHIM to the DXB angular correlation signal in the 
case of a simulated observation with Xenia . In Section [4] 
we measure the spatial two-point correlation function in the 
simulated maps and assess the possibility of characterizing 
the spatial distribution and dynamical state of the WHIM in 
the case of a Xenia observation. In Section [5] we summarize 
and discuss the results of our work. 



2 SIMULATED DATASETS 

In this sec tion we briefly d escrib e the hydrodynamical sim- 
ulation of iBorgani et"al] (|2004|), the WHIM rnodel built 
upon these (jUrsino et al ] I2OIOI : iBranchini et all l2009h and 
the mock X-ray SB maps and spectra obtained from this 
model taking into account the response functions of the 
HARI and ORIS detectors. 



2.1 The Hydrodynamical Model of the WHIM 

The WHIM model considered in this work has been ob- 
tained from the hydrodynamical simulation of lBorgani et al.l 
()2004l ). This simulation assumes a flat {Qm = 0.3, = 0.7) 
ACDM model with a nonvanishing cosmological constant, 
with h = 0.7, Q.b = 0.04, and as = 0.8. The simula- 
tio n was performed using the TREESPH code GAD GET- 
2 (|Springel. White, fc Hernauistll200ll : ISpringelll2005l ) with 
a Plummer equivalent softening, flxed in comoving units, 
epi = 7.5 kpc at 2; < 2, and flxed in physical units 
at 2: > 2. The simulation follows the evolution of 480'' 
dark matter (DM) particles and 480^ baryonic gas par- 
ticles from redshift z = 49 to 2 = 0. The simulation 
box has a side of 192 Mpc. DM and gas particles 

have masses moM = 4.62 x 10® h~^MQ and mgas = 
6.93 X 10* /!,~^M0, respectively. Radiative cooling was per- 
formed under the assumption of optically thin gas with 
pristine composition (76% hydrogen and 24% helium) and 
in coUisional equilibrium. The simulation includes a time- 
dependent photoioni zing UV background th at reionizes the 
Univers e at z ~ 6 (|Haardt fc Madaul [l996l ) . A two- phase 
model (|Sprin gcl fc He rnguistl |2003| ) and a Salpeter initial 
mass function ( Salpc torl 19551 ) were assumed in the simplifled 
star formation model. Galactic winds with typical velocity of 
about 350 km s~^ were assumed to c arry energy, mass, and 
metals produced by stars in the IGM (|Springel fc Hernguistl 
|2003| ). The simulation output consists of 102 boxes, equally 
spaced in logarithmic redshift scale between z — 9 and z — 0. 
This simulation correctly reproduces th e observational prop - 
erties of rich X-ray emitting clusters (|Borgani et al.l [20041 ) 
but overestimates the emission from smaller structures that 
can be identifled with galaxy groups. This overestimate 



probably originates from having assumed primordial com- 
position in the cooling function, instead of updating the 
cool ing rate accordin g to the evolving metallicity of the gas 
(.Bcrtone et al.ll2010l ). However, this effect is not expected 
to affect signiflcantly the emission properties of the WHIM 
gas in regions of lower density, outside virialized structures, 
characterized by an average metallicity well below solar 
IShen et al.ll2010l ) The fact that our WHIM mode l does not 
exceed the DXB constraint (|Roncarelli et al.]|2006l ) corrobo- 
rates this assumption. In fact, the overestimate of the X-ray 
emission from galaxy groups artiflcially reduces the contri- 
bution of the WHIM to the DXB. In this respect, all the- 
oretical prediction on the WHIM detectability presented in 
this work are to be regarded as conservative. 

Metal diffusion is described by particles that actu- 
ally undergo metal enrichment by star formation, and 
not by transfer of metals betwe en neighboring particles 
IWiersma. Schave. fc Smithll2009l ). As a result, the spatial 
distribution of metals may result unrealistically inhomoge- 
neous, in particular for the gas in lower density regions (like 
the WHIM), far away from star forming regions. To circum- 
vent this potential problem, we did not use the metallicity 
calculated self-consistently within the simulation. Instead, 
we assigned a metallicity to all gas particle in each simu- 
lation output by m atching the metallicity- density relation 
in the simulation of ICen fc Ostrikeil (|l999bl ). In practice, at 
a given output, we consider each gas particle with density 
p, and assign a gas metallicity according to the probabil- 
ity function P(Z(p, 2:)) built upon the simulati on outputs of 
ICen fc OstriQ (|l999al ). as fully described in lUrsino et al.l 
( 20101 ). It is reassuring that, with this metallicity prescrip- 
tion, the WHIM model correctly reproduces the observed 
abundances of O vi absorb ers and satisfles the cu rrent con- 
straints for the O VII lines (|Branchini et al.ll2009l ). 

Since we are interested in modeling X-ray observations 
of the WHIM, we focused our analysis on gas at tempera- 
tures higher than 10^ K. We adopted the deflnition of the 
WHIM as the gas with lO'"^ < T < 10^ K and p sC 1000{p), 
distributed along a network of fllaments connecting virial- 
ized structures and in their outskirts. For the purpose of 
characterizing the angular and spatial correlation proper- 
ties of the WHIM, however, we flnd it more convenient to 
consider different classiflcations for the gas phases that we 
describe here for reference. 

In § [21 to study the angular correlation properties of 
the DXB, we flnd it more convenient to divide the gas into 
three phases: 

• The WHIM, which coincides with the previous deflni- 
tion. 

• The non-WHIM. 

• The Total Gas, which represents the sum of the two 
contributions, i.e. all gas with T > 10^ K. 

In § [H to study the spatial correlation properties of the 
line-emitting gas, we consider the following phases: 

• The Total Gas, the gas responsible for detectable line 
emission, irrespective of its thermal state or number density. 

• The Bright-WBIM represents gas regions in which the 
WHIM contributes to the emission of detectable X-ray lines 
(i.e. regions where the WHIM contribution is above the de- 
tectability threshold, described in §[2]3]), regardless whether 
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the non-WHIM contribution is stronger of the WHIM con- 
tribution or not, or whether the non-WHIM contribution is 
above the detectability threshold. 

• The WHIM- dominated gas, represents gas regions in 
which the WHIM contribution to the hne emission is larger 
than that of all other phases and the sum of the WHIM and 
non-WHIM contributions is above the detection threshold, 
even though the WHIM contribution alone is not above the 
detection threshold. 

The Total Gas traces dense and bright regions (usually 
groups and clusters), as well as the less dense regions where 
the WHIM resides (outskirts and bright filaments). The 
Bright-WHTM gas could trace also the brightest regions if, 
model-wise, they are multiphase and have a WHIM contri- 
bution above detectability. The WHIM- dominated gas traces 
only the bright filaments and the outskirts of virialized struc- 
tures, while in the cores of virialized structures the WHIM 
contributes to emission less than the hot and dense gas. 

2.2 Mock X-Ray emission maps 

To study the angular correlation properties of the WHIM 
and assess the possibility of its indirect detection we have 
generated a suite of mock SB maps in the soft X-ray band. 
Here we summarize the procedur e to generate the m aps 
and refer the interested reader to lUrsino et all (|2010l ) for 
a more detailed description. As a first step we consider all 
gas particles within a light-cone, extracted from the hydro- 
dynamical simulation. The light-cone is obtained by stack- 
ing the simulation out put at different epochs out to z — 2 
iRoncarelli et al.ll2006l '). X -ray maps are obtained by consid- 
ering the individual contribution of all the particles within 
the light-cone to the SB. In particular, for every particle 
we obtain the X-ray spectrum from the APEC model of the 
XSPEC0 package to simulate the emission from the plasma. 
Spectra are fully specified by the temperature, density, and 
metallicity of the particle. Here we focus on the energy range 
0.380-0.650 keV, where the WHIM contribution to the DXB 
is thought to be maximum. The lower limit of 0.38 keV al- 
lows us to exclude the strong local C vi line (0.367 keV) 
that would superimpose and dominate the weaker redshifted 
O VII and O viii lines. The upper limit allows us to include 
the O VII triplet and the O viii line at redshift zero, as well 
as redshifted lines of heavier metals like Fe xvii (0.826 keV), 
Ne IX (0.922 keV), and Mg xi (1.352 keV). The use of the 
APEC model is justified under the assumption that the gas 
is in coUisional equilibrium. The goodness of this hypothesis 
is guaranteed by the fact that emission is proportional to 
the square of density. This is also valid for the WHIM con- 
tribution which, therefore, mostly comes from high density 
regions in which the gas is usually hotter and the ionization 
balance is set by collisions only. 

To account for the angular resolution of the detector 
we have binned the particles' contributions to the flux into 
square pixels. The flux in each pixel account for the contri- 
bution from all particles within the pixel itself plus that of 
particles in neighboring pixels weighted according to their 
smoothing kernel. All maps account for the effect of ab- 
sorption due to the neutral hydrogen in our Galaxy. We 

^ http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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Figure 1. Mock Surface brightness in the 380 — 650 eV energy 
band, f.o.v. of 1° X 1° and angular resolution 14 X 14 . Con- 
tributions from the WHIM {top-left) and non-WHIM (top-right) 
phases are shown together with the Total Gas signal (bottom- 
left). Color scale indicates different SB levels in units of pho- 
ton cm~^ s~^. 

modeled absorption according to iMorrison fc McCammonl 
1 19831 ) ■ applying a typical value at high latitude of Nh = 
1.8 X 10^" cm"^ (jMcCammon et al.ll2002l ). Overall we have 
generated 10 simulated maps of the DXB in the energy 
range 380 - 650 eV. We have binned the 1° x 1° f.o.v with 
256 X 256 pixels in order to obtain an angular resolution 
of 14 X 14 comparable to that expected for the HARI in- 
strument. These maps constitute the database that we will 
analyze in §[3l 

Fig. [T] shows a typical mock SB map with a f.o.v. of 
1° X 1° and angular resolution of 14 x 14 . The bottom 
panel shows the flux in the 0.380 — 0.650 keV band con- 
tributed by all gas particles with T > 10^ K (i.e. the Total 
gas). The top-left and top- right panels show the contribution 
of the WHIM and non-WHIM phases, respectively. These 
mock maps also include an isotropic signal cont ributed by 
a a Galactic foreground jMcCammon et al]|2002l 'l and a dif- 
fuse extragalactic background of unresolved AGN, that we 
describe in details in § 13.31 

In Tab. [1] we have listed the mean SB contributed 
by the various components; WHIM, non-WHIM, Galac- 
tic foreground and unresolved AGN. While, as already 
pointed out, the WHIM contribution is consistent with 
observational const raints, the non-WH IM component ex- 
ceeds these limits (jBorgani et al.1 l2004l l. To quantify the 
mismatch is convenient to refer to the 0.65 — 1 keV band 
in which observation a l cons traints have been derived by 
iHickox fc MarkevitchI (|2007l ). The observational value of 
(1.0 ± 0.2) X 10~^^ ergs cm^^ s~^ deg~^ is reassuringly 
larger than the expected WHIM contribution of (3.6±0.3) x 
10~^^ ergs cm~'^ s~^ deg"'^. On the contrary the non-WHIM 
contribution of (8.7 ± 0.5) x 10~^^ ergs cm~^ s~^ deg~^ 
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Table 1 . Contributions to the DXB from the different gas phases 
in the 380 — 650 eV energy band. 



photon cm s ^ deg ^ 


WHIM" 


(1.6 ±0.2) X 10-3 


non-WmW 


(8.9 ± 1.1) X 10-3 


Galactic Foreground'' 


(3.39 ±2.26) X 10-3 


AGN': 


(0.53 ±0.03) X 10-3 



° Mean val ues and standard d e viatio ns calculated from 10 maps. 

Mo del bvlMcCammon et al.l ||2002|). cosmic variance estimated 
from [Gupta et ah I l l2009^JHenlev fc Sheltonl 1 I2OIOI) . 

Double power-la w model (see section[3.3.3l l. errors 
extrapolated from lMoretti et al. 1 1 I2OO3I) . 



largely exceeds the observational constraints as a result of 
the fact that the simulation overestimates the X-ray lumi- 
nosity associated to galaxy groups, as anticipated in the pre- 
vious section. 



2.3 Mock 2D spectra 

To study the spatial correlation properties of the WHIM we 
have simulated data cubes representing a 2D-spectrum sim- 
ilar to that that will be taken by a spectrograph like CRIS. 
For this purpose we have used our WHIM model to obtain a 
"map" consisting of a collection of ID-spectra from contigu- 
ous pixels across a 5.5° x 5.5° area. Here we just summarize 
the main step of the procedure used to g enerate these mock 
spectra and refer the interested reader to lTakei et al.l (|2010l ) 
for a more detailed description. The possibility of performing 
a 3D correlation analysis of the WHIM relies on the pos- 
sibility of detecting its characteristic X-ray emission lines. 
Detection e s timat es using our model have been obtained by 
iTakei et all (j201(]| ) and will be adopted here. 

To produce the spectra we adopt a procedure similar 
to that of the mock SB maps except that here we con- 
sider larger pixels (~ 1.3' x 1.3' corresponding to ~ 50% of 
the CRIS angular resolution) and wider f.o.v. (5.5° x 5.5°, 
roughly corresponding to a mosaic of 25 CRIS pointings). 
The result is a collection of 256 x 256 mock ID-spectra, 
one for each pixel. Furthermore, to account for the energy 
resolution of the instrument, we divide the particle distribu- 
tion along each angular resolution element into radial bins 
of comoving length 3 Mpc h-^ (approximately correspond- 
ing to A£; ~ 1 eV). We end up with 256 x 256 x 448 ~ 
3 x lO'^ mock resolution elements containing spectra con- 
volved with the response matrix of the instrument that 
include the effect of Galactic absorption, modeled as in 
[Morrison fc McCammor] (1 19831), with a c olumn density of 
1.8 x 10^° cm-^ (jMcCammon et al.ll2002l ). 

To reduce the CPU time and memory requirement we 
have considered only gas particles out to z = 0.5, corre- 
sponding to a maximum comoving depth of our correlation 
analysis of ~ 1344 Mpc h-^. In addition, since we are in- 
terested in the O vii and O viii lines only, spectra in each 
bin have been generated using only particles that provide 
contribution in the narrow energy range corresponding to 
that of the two redshifted lines. 

To perform the 3D correlation analysis we consider only 
those resolution elements in which the surface brightness of 



the O VIII line is above 7 x 10~ photon cm- s~ sr- , 
a threshold that guarantees a 5-a d etection in a 1 M s ob- 
servation with the CRIS instrument (|Takei et al.ll2O10[ ). We 
consider the O viii line only, rather than the simultaneous 
detection of O vii and O viii. This choice is justified by the 
fact that ~ 90% of O viii brighter than our selection thresh- 
old is assoc i ated to a detectable O vii line, as shown by 



iTakei et al.l (|2010l ) (Table 1 Model B2, 



1 Ms). The 



visual inspection of the 3D maps corroborates the remark- 
able spatial coincidence between the Q vii+Q viii and the 
O viii-only line-emitting regions (|Takei et al.ll2010l ). 

In the geometry of our problem, the comoving volume 
of the resolution element increases with redshift. Since the 
typical size of line-emitting regions is not expected to evolve 
significantly out to z — 0.5, the risk is of over- or under- 
sampling a line-emitting region, depending on its redshift. 
To circumvent the problem we have grouped the spatially- 
contiguous resolution elements with detectable O viii lines 
into a single object, which we call emitter, and study the 
correlation of such emitters. To assess the robustness of our 
results with respect to the grouping scheme, we have im- 
plemented two different grouping procedures. In the first, 
we link all resolution elements having one side in common. 
In the second, we group the elements that have one edge 
or one vertex in common. We have verified that using ei- 
ther schemes does not significantly modify the results and 
therefore we will only discuss the results obtained with the 
first, more conservative, grouping scheme. The outcome of 
the grouping procedure is a list of O viii emitters character- 
ized by their total O viii line surface (obtained by summing 
over the grouped elements) and spatial position (defined as 
the SB-weighted centre of mass of the emitter). As a final 
remark, we note that the 3 Mpc h-^ binning introduces a 
spurious periodicity along the line of sight which we elimi- 
nate by randomizing the position of the emitter along the 
line of sight, within the resolution element. 



3 ANGULAR CORRELATION ANALYSIS OF 
THE WHIM 

In this section we address the problem of detecting the 
WHIM signature in the angular correlation signal of the 
DXB. We first describe the statistical tool used for the anal- 
ysis, i.e. the two-point correlation function, then we apply 
it to the ideal case of an arbitrarily long exposure time and 
no contribution from AGN and Galactic foreground. Then 
we include all contaminations to simulate a realistic obser- 
vational setup. 



3.1 The angular correlation function 

As shown in Tab. [TJ the Galactic foreground and the non- 
WHIM gas dominate the mean DXB signal in the 380 — 
650 eV energy range. The component we are interested in, 
the WHIM, is largely subdominant (less than 10%) but its 
spatial distribution, from hydro simulations, is markedly dif- 
ferent from that of the other components. Under such cir- 
cumstances, which are reminiscent to that of the CMB sig- 
nal, the best option to infer the presence of the WHIM is by 
analyzing the two-point angular correlation function (ACF) 
of the DXB. 
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We characterize the angular correlation properties of 
the DXB by measuring the angular two-point correlation 
function 'w{9) of the signal in the mock m aps. To measure 
w(6) we use the Landy-Szalay estimator (|Landv fc Szalavl 
1 19931 ) [LS from now on]: 



1 + w{e) 



DD ~ 2DR + RR 
RM ' 



where 6 is the angle between two pixels in the map and 



DD(e) = 
DR{e) = 
RR{e) = 



J2 ^'ki ' 



(1) 

(2) 
(3) 
(4) 



DD represents the number of photon pairs, normalized 
to the total counts, obtained by summing over all photon 
counts gij in two pixels i and j separated by 6. RR is the 
analogous quantity for a map in which the photons are uni- 
formly distributed across the f.o.v (the Random sample). 
DR represents the mixed counts contributed by pixels in 
the real and random dataset. We have verified that the 
LS estimator is equivalent to that used by other authors 
lUrsino fc Galeazzilliooil : iGaleazzi et al.ll2009l '). 

In this work we estimate the ACF of the DXB, the ACFs 
of its components and their cross-terms. For the simple case 
in which the DXB is the sum of two components, A and 
B, with ACF lua and ujb respectively, the following relation 
holds for the ACF of the DXB; 



n^WTiO) = n\wA{0) + n%WB{0) + 2nAnBWAB[d), 



(5) 



where ua and ub are the photon counts contributed by the 
two components, tit = ua -f riB, and ujab is the cross- 
correlation term. 

Uncertainties in the ACF are estimated as the RMS 
scatter among the correlation functions measured among the 
ten, statistically independent, mock maps . 



3.2 Ideal Case 

Are the angular correlation properties of the WHIM signif- 
icantly different from those of the other components? We 
address this question by computing the ACF in the ideal 
case of an arbitrarily long observation with the HARI instru- 
ment, ignoring all additional components like the Galactic 
foreground (which is assumed to be perfectly subtracted) 
and unresolved sources like AGN (that are assumed to pro- 
vide a negligible contribution, as we will show in Section 
[33:3]) . 

For the purpose of detecting the WHIM signature, it is 
useful to compare the ACF of the flux map of the WHIM and 
non- WHIM components. Fig. [2] compares the ACF of the 
WHIM (red, dashed line) with that of the non- WHIM (dot- 
ted, green). Lines represent the mean ACFs of the ten mock 
mapsvvand the band their rms scatter. As expected, the 
non- WHIM provides the dominant contribution to the to- 
tal signal (solid, black line). Even the cross-correlation signal 
(dot-dashed, blue line), estimated by subtracting the WHIM 
and non- WHIM contributions from the total ACF signal, is 
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Figure 2. ACFs of the WHIM {dashed red), non- WHIM [dotted 
green), cross correlation term {dot-dashed blue) and their sum 
(solid black). Amplitudes arc normalized to the total flux. 



higher than that of the WHIM. However the WHIM auto- 
correlation signal is significantly above zero out to S ~ 6'. 
The fact that the WHIM and non- WHIM ACFs have sim- 
ilar shapes is rather disappointing since it means that no 
characteristic feature in the ACF of the WHIM can be used 
to unambiguously infer its presence. On the other hand, the 
small amplitude of the WHIM correlation signal is less of 
a problem, considering that all ACFs in Fig. [2] are normal- 
ized to the total flux. Indeed, the amplitude of the intrinsic 
ACF of the WHIM, i.e. any term w{6) in Eq.O is similar to 
that of the non- WHIM. The difference in amplitude reflects 
the different mean flux of the two phases: a suppression fac- 
tor of ~ 30, with UwHiM and Unon-wHiM the 

y ^TioTi-WHIM y 

WHIM and the non- WHIM SB respectively, is expected be- 
tween the ACF of the WHIM and that of the non- WHIM. 
This is encouraging and suggests that the WHIM correlation 
signal could be extracted from the total one, provided that 
the ACF is estimated in regions where the WHIM emission 
dominates over the non- WHIM one. Such regions are clearly 
seen in Fig.[T]and show that WHIM emission is preferentially 
located at the outskirts of virialized structures whose central 
regions are, on the contrary, dominated by the non- WHIM 
component. 

These mock maps suggest that a possible strategy to 
extract the WHIM autocorrelation signal is to exclude from 
the correlation analysis the pixels associated to the large 
virialized structures in which the SB is usually larger. To 
optimize the cleaning procedure one has to compromise be- 
tween discarding pixels in which the signal is dominated by 
the non- WHIM component, and the need of using as many 
pixels as possible to increase the statistical significance of 
the signal. Our mock maps, in which we can pinpoint the 
contribution of the different components, allow us to cali- 
brate this procedure. We have found that a very simple, yet 
effective, strategy is to exclude pixels with SB above a fixed 
threshold. 

Fig. [3]shows the result of having restricted the estimate 
of the ACF to ~ 50% faintest pixels of the map. The result 
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Figure 3. Same as Fig. |2]after removing the brightest 50% pixels 
from the map. 



demonstrates the goodness of our strategy: the ACF sig- 
nal from the faint regions of the map is indeed dominated 
by the WHIM emission, as expected, and its ACF can be 
regarded as an indirect detection of the WHIM even in ab- 
sence of characteristic signatures in the correlation signal. 
It is worth stressing that the method of discarding bright 
pixels is a very rough one and can be easily improved, for 
example by targeting the pixels associated to known X-ray 
sources. Yet, given the current uncertainties in the WHIM 
model and for the purpose of obtaining an order of mag- 
nitude estimate on the possibility of extracting the WHIM 
correlation signal, we do not try to improve the "cleaning" 
procedure and focus instead on the problem of dealing with 
observational uncertainties. 



pixels. In the second step we assign each pixel the corre- 
sponding (integer) number of photons Nij = int (fijAtexp) 
and we keep track of the difference Pij = fijAtexp — 
Nij G [0, 1). In the third step, we assign the remaining 
Ntot — '^ij photons through Monte Carlo rejection pro- 
cedure that randomly assigns photons to pixels with proba- 
bility P{i,j). Finally, we set the ratio of WHIM/non-WHIM 
photons equal to that of the WHIM/non-WHIM fluxes. 

We did not add Poisson noise to the photon count 
maps since, for long exposure times, the error budget is 
dominated by field-to-field variance rather than shot noise. 
To justify this hypothesis we have evaluated the expected 
Poisson error using the following theoretical expression 
(|Cabre fc Gaztanagal[2009l ): 



AC(r) 



{Nn/NDy+mR/No) 



RR{r) 



(6) 



where Nb. and No represent the number of emitters in the 
real and random catalog. In our case we have used Nr — 
Nd. 

Here we focus on the WHIM signal in the maps from 
which we have removed the 50% brightest pixels. In this 
case we find that field-to-field variance dominates the error 



budget for 



1 Ms. As a result, the ACFs and its er- 



rors computed from the photon maps are identical to those 
computed from the flux maps, in the ideal case (e.g. Fig. [2]). 
However, when texp ~ 100 ks the Poisson noise contribution 
is no more negligible and accounts for ~ 50% of the error 
budget. To account for this effect we have increased the er- 
ror bars on the corresponding maps by a factor \/2. Instead, 
when we consider the Total gas signal in the original maps 
(i.e. with no bright pixel removal) the Poisson noise is largely 
subdominant, even for exposure times as short as 10 ks. 

The same results holds true also when we remove the 
50% brightest pixels from the maps. It is not surprising that 
the ACFs obtained in this case give the same results as in 
the ideal map case (Fig. 



3.3 Realistic Case 

We now repeat the analysis of § 13.21 using more realistic 
mock maps which account for instrumental effects (expo- 
sure time and collecting area) and for the presence of quasi- 
isotropic fore/backgrounds (Galactic foreground and unre- 
solved AGN). 



3.3.1 Discrete photon counts 

In a realistic observational setup with finite exposure time 
and collecting area, one deals with photon counts. In order 
to transform the maps of § 13.21 into photon counts maps, 
we multiply the flux by the exposure time and the effective 
area of the detector. We consider observations with t^^p in 
the range [0.1-1] Ms and a detector of effective area A = 
10^ cm^, to mimic a long/very long exposure with the HARI 
ccd. 

To generate the discrete photon maps we adopt a 
four-step procedure. As a first step we compute the to- 
tal (integer) number of photons in the map as Ntot = 

int ^X^ij fijAtexp^ , where fij is the expected flux in the 

pixel with coordinates (i, j) and the sum runs over all the 



3.3.2 Galactic Foreground 

Realistic DXB maps must account for the local contribu- 
tions to the signal, i.e. for the Galactic foreground, which 
we model with two components: the Local Bubble and the 
Galactic Halo. The Local Bubble emission, generated by the 
hot gas within ~ 200 pc, is conveniently modeled as an un- 
absorbed thermal component with T ~ 10^ K and solar 



meta l abundances |S nowd cn ct al. 1992 



I2OO2I : ICaleazzi et aljr2007l : ICupta et al.ll2009l ). In our maps 



iMcCammon et aP 



it is included using the APEC model in the XSPEC pack- 
age. The Galactic Halo, generated by thermal emission of 
gas within our Galaxy, is well described by an absorbed 
thermal component with T ^ 2 x 10^ K JSnowden et al.l 
1 19971 : iMcCammon et all |2002| : iGaleazzi et all bOOTl ). It is 
added to our maps by combining an APEC and a WABS 
model in XS PEC. Both models use th e parameters listed 
in Tab. 3 of IMcCammon et all (|2002l '). The resulting SB 
of the Galactic foreground in the 380 — 650 eV band 
is listed in Tab. [T] This galactic foreground model is 
effective since it matches the characteristics of the ob- 
served flux attributed to Galactic sources. Yet it is not 
accurate since it ignores the contribution from the Solar 
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Wind Charge Exchange (SWCX) whic h has been shown 
to ou t shine t hat of the Local Bubble ( Koutroumpa et alj 



2006^ 2007; 'Hcnlcv & Shclton' '2008'; 'Gupt a et all l2009l : 
Kout roumpa ct al. 2009.: .Hcnlcv fc Sholtoni,2010l ). The vari- 
able nature of the SWCX and the lack of strong observa- 
tional constraints makes it difficult to model its contribution 
to the DXB. Fortunately, the very local origin of the SWCX 
excludes the presence of small (i.e < 10') angular fluctua- 
tions so that its contribution to the Galactic foreground can 
be effecti vely modeled with an isotropic component which, 
following [Gupta et all (|2009h . we obtain by boosting up the 
contribution of the Local Bubble and Galactic Halo. 

Since angular variations of the Galactic foreground are 
expected to be significant only on angular scales much larger 
that those relevant for our correlation study, we model it as a 
purely isotropic component that contributes a Nf photons 



in each pixel of the maps. For te 



1 Ms, Nf = 51.76. 



Therefore, we obtain the number of Galactic foreground 
photons in the maps by Monte Carlo sampling a Poisson 
distribution with mean Nf- 

The addition of the Galactic foreground significantly re- 
duces the possibility of detecting the WHIM signature in the 
DXB in two ways. First, it decreases the amplitude of the 
ACF signal. According to Eq. [5] the expected decrement is 



50% in the case of te 



1 Ms. The effect on the signal- 



to-noise error, however is less dramatic since the presence 
of a bright isotropic foreground greatly reduces the field-to- 
field scatter which dominates the noise. In fact, the overall 
reduction of the statistical significance of the correlation sig- 
nal is just ~ 20%. Second, and more serious, is the impact of 
the Galactic foreground on our capability of extracting the 
WHIM signal from the maps. The addition of a bright fore- 
ground outshines the correlation signal in the faint pixels. 
As a consequence, removing the brightest 50% pixels results 
in shot-noise dominated maps with zero ACF at all angles, 
as we have verified. 

A simple solution to this problem is that of removing 
a model Galactic foreground before analyzing the ACF of 
the maps. Foreground removal would be then performed by 
simply subtracting the Galactic foreground model (a ran- 
dom counts map characterized by its mean intensity) to the 
photon count map. Since we assume that random counts 
are Poisson distributed, the subtraction is performed as fol- 
lows. First we Monte Carlo generate the number of random 
counts in the pixel, Fij . Then we compute the residual map 



Ri 



OifT„ < Fu 



Since this procedure removes less photons from the total 
map than those in the foreground map Fij, we randomly 
remove the remaining photons from the pixels with Rij > 0. 
Finally, we remove the 50% brightest pixels from this resid- 
ual map. The result is displayed in Fig. 3] The panels show 
the same f.o.v. as Fig. [1] The extended, black areas corre- 
spond to the regions in which the 50% brightest pixels have 
been removed from the map. Isolated black pixels in the 
low-count regions are those with Rij =0. 

This procedure to remove the foreground is very sim- 
plistic but, as shown in the top panel of Fig. [5] it is good 
enough to extract the WHIM autocorrelation signal from a 
tsxp = 1 Ms observation, since the ACF below S ~ 3' is 
significantly above zero and dominated by the WHIM com- 
ponent (red dashed curve). 

As in the previous section, we do not add Poisson noise 



to the WHIM and non- WHIM counts. This choice is justified 
because, in addition to the fact that cosmic variance domi- 
nates over the Poisson noise of the signal for f te^p = 1 Ms, 
as discussed previously, the Poisson noise is now dominated 
by the Foreground, not by the signal, which is included in 
our maps. 

Given the uncertainties in the WHIM model adopted, 
these results are to be regarded as indicative. Yet, they serve 
to illustrate an important point: the capability of extract- 
ing the WHIM autocorrelation signal from the DXB maps is 
mainly hampered by the presence of a bright Galactic fore- 
ground. With this respect, our problem is analogous to that 
of extracting a weak signal in presence of a strong random 
noise. In this case, the signal-to-noise can be improved in 
two ways. 

Since, in the large noise limit, signal-to-noise ratio scales 
as texp, the simplest option is to increase the exposure time. 
The effect of decreasing texp by a factor 10 is illustrated 
in the bottom panel of Fig. [S] In this case the ACF signal 
has been washed out everywhere but at sub-arcmin angular 
separations. Varying the exposure time should not affect the 
ACF signal but only its uncertainties, if they are dominated 
by Poisson noise. However, the two panels of Fig. [5] show 
that this is not the case. We have verified that this apparent 
inconsistency reflects a systematic bias introduced by our 
procedure to remove the foreground from a map in which 
the signal is dominated by the foreground itself. The shorter 
the exposure time, the less efficient the foreground removal, 
the more Poisson noise dominates the residual map. The 
net results is to wash out the ACF signal in the maps with 
shorter exposure times. 

To assess the accuracy of our simple strategy to remove 
the foreground we have considered an additional procedure 
in which we remove a constant signal with amplitude equal 
to the mean foreground, but with no Poisson noise included, 
with the caution of avoiding negative photon counts. The 
results are basically the same as with the original procedure. 

Rather then increasing the exposure time, a better al- 
ternative is to improve the map cleaning efficiency. First 
of all, one can perform a more effective exclusion of the 
non- WHIM signal by excluding (or filtering out) the signal 
in correspondence of known X-ray sources. Second of all, a 
better extraction of the WHIM signal can b e obtained by 
convolving the m ap with Wiener- like filters (|Wieneilll949l : 
IPress et al]|l992l ). Of course Wiener filtering relies on a- 
priori modeling of the WHIM ACF and Galactic foreground 
noise. However, possible inaccurate values of the Wiener Fil- 
ter resulting from mo del uncertainties would introduce sec- 
ond order errors only (|Rvbickv fc Presslll992l V 

Our analysis, therefore, clearly shows that long expo- 
sures are required to extract the WHIM autocorrelation sig- 
nal from the ACF of the DXB. However, our considerations 
suggest that the value of texp ~ 1 Ms is probably too gener- 
ous and that a significant detection might be possible with 
shorter exposu res. This conclusion is corroborated by the re- 
cent results of ICaleazzi et all (|2009l ), who studied the ACF 
of the DXB in a set of XMM-fields with texp ranging from 
50 ks to 600 ks. After removing the contribution from known 
sources and making no attempt to subtract the Galactic 
foreground, they d etect a positive auto correlation signal in 
all fields (Fig. 2 of ICaleazzi et al1l2009l ) that they interpret 
as the signature of the WHIM. 
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Figure 4. Residual maps of WHIM (top-left), non-WHIM {top- 
right), and Total gas (bottom-left) in the 380 — 650 eV energy 
band for a f.o.v. of 1 X 1 deg^, effective area of 1000 cm'^, and 
texp = 1 Ms. The color scale is in photon-number units. 



Finally, we comment on the error contribution from 
the expected instrumental noise. The baseline of Xenia 
CCD detector is similar to that of Suzaku and relies on 
being placed on a low- Earth orbit. The rational behind 
this choice is to guarantee a more stable background than 
that in highly elliptical orbits where Chandra and XMM- 
Newton have been placed. The Suzaku instrumental back- 
ground is 

FI and ~ 5 x 10~^ counts kev" 
BI, i.e. 4 ti mes and 9 times s maller than the cosmic X-ray 
background (|Tawa et al By comparison instrumental 

backgrounds for XMM-MOS and for XMM-PN are 2 times 
and 5 times smaller than the cosmic X-ray background. 
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Figure 5. Comparison between the ACF of total map {solid 
black), WHIM component {dashed red), raon-WHIM component 
{dotted green), and the cross-correlation {dot-dashed blue). All 
ACFs are normalized to the total flux after removing the bright- 
est 50% of the pixels. The plots refer to the case of a texp = 1 Ms 
observation (top) and a texp = 100 ks observation {bottom). 



3.3.3 Resolved and unresolved AGN 

We also include the contribution of resolved and unresolved 
AGN in the mock maps. 

Resolved AGN can be easily eliminated from the maps 
by excluding the contaminated pixels. However, given the 
limited number of expected objects, the effect is rather 
small. For texp ~ 1 Ms, Xema is expected to resolve ~ 2000 
AGN per square degree above the detection thres hold Sxenia 
of 10 "^'' erg cm~^ s"^ in the 0.5 - 2 keV band (|Piro et al.1 
I2OO9I ). Under the somewhat simplistic assumption that each 
resolved AGN occupies just one pixel, the fraction of ex- 
cluded pixels is just 3% of the total, much smaller than the 
50% fraction of bright pixels that were excluded to remove 
the non- WHIM contribution. It is no surprise that remov- 
ing this extra 3% of pixels has no appreciable impact on 
our result. It is also clear that a more sophisticated model, 
in which one account for the instrumental PSF, would not 
affect our results either. In fact the situation is even more 



favorable, since most of the AGN will be associated to mas- 
sive structures contributing to the non- WHIM signal corre- 
sponding to the already removed 50% brightest pixels. 

To estimate the effect of unresolved object, we have in- 
cluded their contribution to the DXB by generating a pop- 
ulation of unresolved sources matching the number counts 
and redshift distribution of faint AGN observed by Chan- 
dra and using N-body simulation to model their spatial cor- 
relation properties. More specifically, we generate a mock 
population of AGN as follows. 

First we assume that mock AGN have the following 
number counts: 

Tir/ c\ J ) ^ SLim ^ 5* *C Schandra 

I , LX. uChandra ^ J ^ OXema 

where Schandra = 3 X 10~^^ erg cm~^ s"'^ is similar 
to the detection limit in the Chandra Deep Fields and 
SLim ~ 6.5 X 10"^** erg cm~^ s~^ is a lower limit so that 
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the function does not diverge at infinity, fn practice we 
require tfiat mock AGN counts matcfi tfie Chandra ones 
above Shim (|Moretti et al.|[2003h and assume tiiat tfiey fol- 
low the usual Euclidean slope below SLim- The reason of 
assuming a Euclidean slope rather than extrapolating the 
one of the resolved AGN is to maximize the contribution 
of unresolved objects to the DXB. To assign flux to the 
mock unresolved AGN we sampled the probability distri- 
bution obtained from Eq. [7] between Sxenia to SLim- The 
total number of object is determined by requiring that un- 
resolved AGN contribute to ~ 12% of the Cosmic X-Ray 
Background (assuming a Cosmic X- Ray Background flu x 
of 7.5 X 10"^^ erg cm"^ s"^ deg"^ iMoretti et al]|2003h . 
In Tab. [l] we report the flux of the unresolved AGN in 
the 0.38 — 0.65 keV energy band assuming a spectral slope 
r = 1.4. We end up by generating ~ 5 x 10** objects per 
square degree with SLim ^ S < Sxenia- 

To determine the redshift distribution of these 

dz 

objects we in t egrat e the AGN luminosity function of 
iHopk ins et al] ()2006l ) specified in three different redshift 
bins, out to z = 3. 

Finally, we determine the spatial location of the mock 
AGN from the outputs of a coUisionless N-body experiment 
under the assumption that AGN are located at the centre 
of dark matter halos. In practice, we extract 10 independent 
light-cones with f. o.v. of 1° x 1° from the publicly available 
GIF simulations (| Jenkins et al.|[200lh . Light cones are ob- 
tained by stacking all available outputs of the 479 Mpc h^'^ 
computational cube out to z = 3. Within each cone we adopt 
a Monte Carlo rejection procedure to select a subsample of 
5 X 10* dark matter halos with mass > 1.4 x 10^^ Mq h~^, 
enforcing the required X-ray fluxes are randomly as- 
signed according to Eq. [7] and the angular position of the 
AGN within each light-cone is obtained by projecting the 
spatial position of each selected halo on the 1° x 1° f.o.v. 

By construction, mock AGN have the same autocor- 
relation function of the dark matter halos, but no cross- 
correlation with the diffuse gas, since gas particles and dark 
matter halos have been taken from different simulations. 

This inconsistency, however, is hardly an issue since, as 
we have verified, the angular autocorrelation function of the 
unresolved AGN is consistent with zero within the errors. 
We conclude that refining the AGN model to account for 
the spatial correlation with the gas distribution would not 
change our results. Unresolved AGN fainter than SLim pro- 
vide no significant correlation signal and can be regarded 
as an additional, wea k isotropic signal t o be added to the 
Galactic Foreground (|Plionis et al.ll2008l '). 

Another potential class of sources that may contribute 
to the DXB is that of normal galaxies, thanks to their X-ray 
coronae and X-ray binary population. However, at the Xenia 
detection threshold for 1 Ms exposure, 10~^^ erg cm^'^ s^^, 
their contributio n is largely dominated by tha t of the AGN, 
as shown e.g. bv lHickox fc MarkevitchI (|2006l ). To indepen- 
dently check the validity of this conclusion we have extrap- 
olated the Log iV/Log S of normal galaxies computed by 
iGeorgakakis et all (|2006l ) (Fig. 3 of their paper) down to 
6.5 X 10~^® — 1 X 10"^'' ergs s~^ cm~^ and computed their 
contribution to the diffuse background of Xenia . In this 
exercise we had to transform the original flux in the [0.5 — 
2] keV band into the X-ray band considered in this paper, 
[0.380 — 0.650] keV. For this purpose we have assumed that 



the typical 0.380 — 0.650/0.5 — 2 keV flux ratio for a normal 
galaxy is identical to that of the Milky Way. A choice that 
probably provides an upper limit since the the many Galac- 
tic lines that contribute to the flux in the [0.380 — 0.650] keV 
are redshifted out of this range in the spectrum of a distant 
object. We flnd that the contribution of unresolved galaxies 
is 2.78 X 10"* phot s"^ cm ^ sr ^, i.e. a factor of ~ 2 smaller 
than that of the AGN. Had we assumed a more conserva- 
tive, single powe r law A'^(> S) cc S~^' ^ for t he unresolved 
AGN, as e.g. in iHickox fc MarkevitchI (|2006l '). the relative 
contribution of normal galaxies would have been larger. 

At lower fluxes the contribution of normal galaxies is 
expected to dominate ove r the AGN, but not over th at of 
the WHIM, as discussed in lHickox fc MarkevitchI (|2007^ and 
detailed in Table [T] 



4 TRACING THE SPATIAL DISTRIBUTION 
OF THE WHIM 

In this section we explore the possibility of performing 
2D spectroscopy to characterize the 3D distribution of the 
WHIM. For this purpose we take the CRIS instrument as 
reference and use our mock 2D spectra to assess how well one 
can trace the 3D distribution of the WHIM using its char- 
acteristic O VII and O viii emission lines. The angular and 
energy resolution of the instrument limits the ability of trac- 
ing the angular and redshift distribution of the WHIM, re- 
spectively. The exposure time determines the sampling den- 
sity of the line-emitting gas. Here we consider the case of 

texp = 1 Ms. 

The number of WHIM detections expect ed with an 
instru ment like CRIS has been estimated by iTakei et al.l 
(|2010h . Their analysis shows that the minimum line sur- 
face brightness required for a 5a detection is 0.07 pho- 
ton cm~^ s~^ sr~^and that the expected number of simulta- 
neous O VII and O viii detections is 639 per square degree, 
426 of which are contributed by the WHIM (see Tab. 2 of 
iTakei et al.ll2010l ). In fact, the actual number of detections 
is expected to be significantly larger, since these estimates 
were obtained considering gas with z < 0.5, whereas the 
WHIM mass fraction is still large up to 2: ~ 1. 

These detection estimates, however, have been obtained 
for an angular resolution (2.6' x 2.6') larger than that of our 
mock 2D spectra (1.3' x 1.3'). As a consequence, we should 
revise detection estimates to account for the larger number 
of angular resolution elements available and for the larger 
threshold (0.22 photon cm~^ s"'^ sr~^) required for a 5a 
detection. However, our goal is to detect line emitters that 
were obtained from a grouping procedure. It turns out that 
typical angular size of the emitters is closer to 2.6' x 2.6' 
th an to 1.3' X 1 .3' an d we decide to keep the same threshold 
as iTakei et al.1 (|2010l ') to select detectable emitters. The fact 
that the number of d e tecta ble emitters is close to the esti- 
mates of iTakei et al.l (|2010l ') confirms the goodness of this 
choice. 

To characterize the spatial distribution of the line emit- 
ters we use the spatial two-point correlation function (,{r), 
which we compute using the same estimator used for the 
angular correlation function (see Section I3.1|l , generalized 
to the three dimensional case. The analysis is performed 
both in real and redshift space, i.e. by considering both co- 
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moving distances and "observed" redshifts. Finally, as al- 
ready pointed out, since for our emitters the presence of a 
detectable O viii line often guarantees the simultaneous de- 
tection of the O VII line, in our analysis we will only consider 
O VIII emitters. 



4.1 Redshift distribution of line emitters 

Before studying the spatial correlation properties of the 
O VIII emitters, let us analyze their redshift distribution out 
to the comoving distance dc — 1300 Mpc h~^, corresponding 
to the maximum redshift z = 0.5. In this redshift range the 
WHIM mass fraction is expected to increase by only ~ 10%, 
whereas the mass fra ction of the hot gas {T > 10^ K) should 

If O VIII lines 



mcrease by ~ 50% (|Cen fc Ostrike: 



gasJT ; 
3[200i). 



were mostly produced by hot gas in virialized structures we 
should observe a significant evolution in the number density 
of O VIII emitters. On the contrary, if emission were domi- 
nated by the WHIM, the emitter number density should re- 
main constant. Here we briefly remind the reader of the def- 
inition adopted in Section [2. II The Total Gas consists of all 
gas with temperature above lO'' K. The Bright-W}11M rep- 
resents where the WHIM contributes to the X-ray emission 
is detectable (regardless of the non-WHIM). The WHIM- 
dominated gas consists of gas regions in which the WHIM 
contribution to the line emission is larger than that of all 
other phases. 

Fig. [6] shows the cumulative density distribution of the 
O viii-line regions in the simulated light-cone as a function of 
the comoving distance dc- Once we ignore the local region 
{dc < 300 Mpc h~^) dominated by a few structures and 
characterized by a large variance, we see that the comoving 
number density of all O viii emitters (green-dotted curve) 
is not too different from that of a volume limited sample, 
A''(< dc) oc d'l (blue dash-dotted). Deviations from this law 
exist for dc > 500 Mpc h^^. However, the fact that the slope 
of the Bright-WWIM curve (black, solid) is similar to that of 
the Total Gas suggests that they are not due to evolutionary 
effects. In fact, the similar shapes of the curves indicates that 
in both cases the emission is dominated by diffuse material, 
whose mass fraction hardly changes in this redshift interval, 
rather than by virialized gas. This result simply confirms 
that, due to the unfavorable ionization balance, the bulk of 
the O viii-line emitters does not trace the hot {T > 10^ K) 
gas associated to galaxy clusters. Instead, they trace the 
gas located in the outskirts of virialized objects. Deviations 
from the A''(< dc) oc dc curve simply derives from the fact 
that the volume of the resolution elements increases with the 
distance. This means that two nearby line-emitting regions 
along the same line of sight would be counted as two different 
emitters when they are close to the observer but would count 
as a single object at large distances. 

The very fact that the cumulative number density of the 
WYIIM- dominated (red, dashed) emitters decreases faster 
with the distance than that of the Brig/ii- WHIM emitters 
corroborates this hypothesis. To understand this point, one 
has to keep in mind that the Brig/ii- WHIM gas is found in 
the same regions, the outskirts of virialized objects, as the 
"WBIM- dominated. The main difference is that the latter is, 
on average, considerably fainter, i.e. its lines are closer to 
the detection threshold. Therefore an extended region that 
looks like a single spot, if regarded as Brig/ii-WHIM, would 
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Figure 6. Cumulative distribution of O VIII emitters normalized 
to that of a population of objects with constant number den- 
sity oc d^. Bright WHIM (solid, black), WHIM- dominated (red, 
dashed) and Total gas (green, dotted) emitters. The blue, dash- 
dotted line is a reference distribution for a volume limited sample 
(A'(< dc)ocdl). 



have a more patchy appearance when regarded as WHIM- 
dominated gas. The grouping procedure might then extract 
several WYLIM- dominated emitters from a region in which 
only a single Bright-WBIM emitter is found. The magnitude 
of the oversampling decreases with the redshift as demon- 
strated by the fact that in Fig. [6]the total number of Bright- 
WHIM emitters matches that of the WHIM- rfommaied emit- 
ters at dc > 1000 Mpc h"\ 

To summarize, the cumulative redshift distribution of 
the O VIII emitters suggests little or no evolution in the 
range z G [0,0.5]. Deviations from a population of objects 
with constant number density and discrepancies between 
WYIIM- dominated and Bright-WHIM emitters can be un- 
derstood in terms of geometry effects. As fully explained in 
the following section, the errors shown in Fig. [6] are not es- 
timated as scatter between different samples since we have 
only one catalog, instead we assumed a Poissonian distribu- 
tion of emitters, and the errors are simply the square root 
of emitters' number. 



4.2 Two-point correlation function in real space 

Fig. [7] shows the spatial two-point correlation function, ^(r), 
of the detectable O viii-line emitters as a function of the 
comoving pair separation. The correlation function is mea- 
sured in real space, ignoring peculiar velocities. The curves 
refer to the three families of emitters: Bright WHIM (black, 
solid), WBIM- dominated (red, dashed) and Total gas (green, 
dotted). The correlation functions are computed in comov- 
ing bins of 3 Mpc h~^ to account for the instrument energy 
resolution. Ideally, we would like to estimate errors from 
the scatter among the mocks, as in § [S] However, in this 
case we only have one mock catalog at our disposal. For- 
tunately, given the large volume covered by our mock 2D 
spectra and the limited number of detectable emitters, we 
expect that the errors in the ^(r) are dominated by sparse 
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sampling rather than by field-to-field variance. To model 
shot noise we use Eq.[6l as we did in Section [3. 3. II this time 
adopting Nr/Nd = 3. In this work we rely on Eq. [6] rather 
than computing errors using more sophisticated techniques 
(bootstrap resampling, jackknife) since the main goal here 
is to obtain an order-of-magnitude estimate for the errors to 
Eissess our fiiture capability of characterizing the correlation 
properties of the WHIM. In practice we will only trust dif- 
ferences that are significantly larger than the l-o width of 
the error-bands in Fig. [7] 

The main result of our analysis, shown in Fig.[71 is that 
O VIII emitters, whether or not contributed by the WHIM, 
are expected to be spatially correlated and their two point 
correlation function will be detected with high statistical 
significance by next generation experiments like Xema. The 
large signal-to-noise ratio in the expected £, (r) demonstrates 
the robustness of this prediction which cannot be too much 
affected by uncertainties in the WHIM model and approxi- 
mations in our error estimates. 

All two-point correlation functions in Fig. [7] are well 
approximated by a power law with slope 7 = —1.7 ± 0.1 
in the range r = [4.5,12] Mpc h~^, steepening at larger 
separations. Such slope, determined via minimization, is 
common to many cosmic structures, from galaxies to clusters 
and simply shows that both the O viii emitters and virialized 
objects trace the underlying distribution of Dark Matter, as 
expected in a standard CDM scenario. 

In particular, the correlation length of all detectable 
emitters (the Total gas) ro = 4.0 ± 0.1 Mpc h~^ is consis- 
tent with that of opticall y selected galaxies as measured in 
sever al redshift catalogs jZehavi et al.|[20o3 : iHawkins et al.l 
I2OO3I ). Such coincidence does not necessarily indicate a 1:1 
relation between galaxies and O viii emitters. Instead, it 
shows that both types of objects trace the same dark matter 
structures at separations > 3 Mpc h~^, significantly larger 
than the virial radius of galaxy groups and clusters. 

One interesting feature of Fig. [7] is that the correlation 
of the emitters contributed by the WHIM {Bright-WHIM 
and 'W}ilM- dominated) is stronger than that of the Total 
gas, whose line emission is contributed by gas in groups 
and clusters. This result is counter-intuitive since hydro- 
simulations show that the WHIM preferentially traces the 
filamentary structure of the cosmic web whose spatial cor- 
relation is weaker than that of virialized objects. This re- 
sults can be understood as the combination of two effects: 
a selection effect, derived from having impos ed a given SB 
threshold, and the peak-biasing phenomenon (|Kaiseilll984h . 
which results from sampling the peaks of the underlying 
mass density field. In the case we are considering, the O viii- 
line associated to the WHIM in filamentary structures is too 
weak to be detected. Only the WHIM in the outskirts of 
large, virialized structures is O viii-bright enough to be de- 
tected. In other words: with our selection we only sample 
Brig/it- WHIM and W}HM- dominated emitters associated to 
galaxy clusters, whereas Total gas emitters are preferentially 
associated to the external regions of smaller, galaxy group- 
size structures. Galaxy clusters are associated to the highest 
peaks of the mass density fields, higher than those associated 
to galaxy groups. The high er the peaks, the larger their two- 
point correlation function (|Kaiseill 19841 ) ■ which explains the 
stronger correlation of WHIM-contributed O viii emitters. 

We stress that the explanation above, though generally 



valid, has to be regarded as semi-quantitative since O viii 
emitters are not directly associated to virialized structures, 
but to their outskirts, i.e. regions beyond their virial radius. 
For this reason it is not possible to interpret our result in 
terms of a halo biasing model, and infer the typical prop- 
erties of dark matter halos from the correlation function of 
the O VIII emitters. 

In Fig. [7] we also notice that the WHIM- rfommaied emit- 
ters are more correlated than the Bright-W}11M ones. The 
reason for this difference is more subtle. We have seen that 
both types of emitters sample the same regions (the outskirts 
of galaxy clusters) and so they should have the same correla- 
tion function. However, as we have discussed in Section [4. 11 
the W}iTM- dominated emitters oversample these regions due 
to their patchy appearance. The result of oversampling the 
very same peaks of the density fields is to increase the cor- 
relation of the WBIM- dominated with respect to that of the 
Bright-WniM emitters. 

The small scale (^5 3 Mpc h~^) fiattening of the cor- 
relation function is an artifact. It refiects the fact that the 
energy resolution of the instrument effectively acts as an 
anisotropic smoothing filter with radial scale of 3 Mpc h~^ 
which erases correlation on the corresponding scales. Once 
again, the different behavior of the WYilM- dominated emit- 
ters reflects their oversampling of the density peaks, which 
is particularly effective at small separations. 

Fig. [7] shows the correlation function measured in the 
full light-cone. Since the mass fractions of different gas 
phases evolve differently we should check whether their cor- 
relation functions evolve in a different way too. To inves- 
tigate this issue, we have measured the spatial correlation 
function of the O viii emitters slicing the interval z — [0, 0.5] 
into six bins. The results are shown in Fig. [H] 

There is little evolution in the correlation properties 
of the emitters regions in the interval explored, with the 
exception of the nearest and outermost slices. At z ~ 0.5 the 
WYIIM- dominated regions appear to be more correlated on 
larger scales than at lower redshifts. The simple explanation 
is, again, related to selection effect and peak biasing. At 
large redshift one detects stronger O viii emitters associated 
to higher density peaks and the amplitude of the correlation 
function increases accordingly. On the opposite, at 2; ~ all 
types of emitters are less correlated than at higher redshifts. 
The effect is particularly visible at large separations but its 
statistical significance is not large. The point is that the 
volume associated to the first slice of the cone is rather small 
and does not allow to fairly sample scales 10 Mpc . 
Therefore, the error budget is likely to be dominated by 
cosmic variance rather than shot-noise errors. Adding the 
cosmic variance contribution would increase the size of the 
error-bands in Fig. [7| and reduce the statistical significance 
of the effect. 

4.3 Correlation analysis in redshift space 

From the analysis of the X-ray spectra one infers the 3D 
distribution of the O viii emitters in redshift space, not in 
real space. The presence of peculiar velocity spoils the 1:1 
relation between redshifts and comoving distances and in- 
troduces systematic distortions in the spatial clustering of 
the emitters and in their correlation properties. The net re- 
sult is that the two-point correlation function measured in 



Expected properties of the Two-Point Autocorrelation Function of the IGM 13 



10 F 



O.I - 



0.01 




I 10 
r (Mpc h"') 



30 



Figure 7. Real-space correlation function for ftrig/itWHIM (solid 
black), WHIM- dominated (dashed red), and Total gas (dotted 
green) O VIII emitters. Shaded areas represent Poissonian un- 
certainties. 
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Figure 8. Real-space correlation function for bright-WKIM, 
(solid black), WHIM- dominated (dashed red), and Total gas (dot- 
ted green) O VIII emitters in six redshift slices. The redshift range 
is indicated in each panel. 



redshift space is systematically different from the true one, 
measured in real space. The mismatch can be used to infer 
the properties of the peculiar velocity field, i.e. to character- 
ize the dynamical state of the emitters. 

Peculiar motions can be roughly divided in two cat- 
egories. Incoherent motions, which typically dominate at 
small separations, and represent random velocities within 
virialized structures. They cause these structures to appear 
elongated along the line of sight in a characteristic type of 
distortion commonly known as "Fingers of God". Coher- 



ent motions characterize large scale structures growing in 
the linear or quasi-linear regime and typically represent in- 
fall motions into mass overdense regions (or outflow from 
underdense regions). They apparently increase the density 
contrast associated to these regions along the line of sight. 

For a collisionless fluid, like the galaxies in the universe, 
measuring these motions allows to infer important proper- 
ties of the underlying mass density held: the amplitude of 
virial motions allows to infer the typical mass of virialized 
structures in the sample, while the size of coherent motions 
allows to constr ain the g rowth rate of density fluctuations 
(,Hamihon.l993 ;[G uzzo et al. 2008,) fro m which one can mea - 
sure the mass density parameter Qm (iPeacock et al.l 1200 ll ) 
or test alternative gravity models (|Lindeijl2005l ). For the case 
of a collisional dissipative fluid, like the O viii emitting gas, 
one cannot derive these parameters from redshift distortions. 
However, it is still possible to infer the dynamical status of 
the gas which still brings some important information on 
the typical environment in which the WHIM is found. A 
strong Finger of God signature would indicate that the gas 
is preferentially found within virialized environments, while 
a dominance of the compression effect by coherent motions 
would indicate that the WHIM is preferentially found in 
regions of moderate overdensity moving toward virialized 
regions or outflowing from low density environments. 

Peculiar velocities cause distortions only along the line 
of sight, and not on the transverse direction. One can there- 
fore quantify their amplitude by comparing the clustering 
of objects in the directions parallel and perpendicular to 
the line of sight. For this purpose it is convenient to com- 
pute the two point correlation function after decomposing 
pair separations r into the line of sight component, n, and 
perpendicular to it, Vp. The resulting correlation function 
^(rp,7r) is readily computed using the usual LS estimator. 
The three panels in Fig. [9] show the two-point correlation 
function ^{rp,n) of the three types of O viii emitters, pro- 
jected onto the Vp, n plane. Different colors indicate different 
iso-correlation contours, whose amplitude is indicated in the 
color scale. 

With no peculiar velocities, clustering properties are ex- 
pected to be isotropic and iso-correlation contours would 
be circularly symmetric. Deviations from the circular sym- 
metry would then indicate the presence of peculiar-velocity 
induced anisotropics in the spatial clustering of the emit- 
ters. The stronger the deviations, the larger the amplitude 
of peculiar velocities. As anticipated these distortions can 
roughly be divided in two categories. Incoherent velocities 
associated to Fingers of God type distortions cause the iso- 
correlation contours to be elongated along the tt direction. 
Being associated to small-scale motions they constitute the 
dominant type of distortions at small values of rp. Coherent 
motions associated to enhancements of density contrasts in- 
duce a compression of the iso-correlation contours along tt. 
Since they are associated to large-scale motions, they dom- 
inate the distortion pattern at large values of rp. 

The iso-correlation contours shown in Fig. [9] display a 
non- negligible compression for Vp ^ 6 Mpc h~^ whose mag- 
nitude does not seem to depend on the type of emitter. This 
result suggests that on large enough scales, O viii emitters 
trace the same underlying velocity held, characterized by 
linear, coherent flows, irrespective of their WHIM content. 
At smaller scales the iso-correlation contours appear to be 



14 E. Ursino et al. 



elongated along the line of sight which characterize Fingers 
of God- type distortions. In fact, these distortions appear 
to be more prominent for WYiliA- dominated emitters, while 
they seem to be absent in the Total gas case. 

To make our analysis more quantitative and assess 
the statistical significance of these distortions, we have fit- 
ted the measured (,{rp,Ti) with an analytic model that ac- 
counts for both large-scale coherent motions and small- 
scales by expanding ^(r;3,7r ) in Legendr e polynomials out 
to the quadrupole term (jKaiseil 1 19871 : lHamilto"nl 1 19921 : 
ICabre fc Gaztafiagal [ioool ). This model is fully character- 
ized by two parameters: an "elongation" parameter g\2, 
which quantifies the strength of Fingers of God-type dis- 
tortions, and a "compression" parameter, /3, which quanti- 
fies the strength of linear-type distortions. For a coUisionless 
fluid, a"i2 measures the typical strength of incoherent mo- 
tions while /3 quantifies the amplitude of coherent motions, 
which is proportional to the growth of density fiuctuations 
/: /? oc /(r2„i) oc fim^*, where is the mean mass density 
in units of critical density. For O viii emitters such relations 
are not valid anymore. However, one can still measure ai2 
and /3, and compare them with the null case to detect the 
presence of either type of motions, or to compare the rela- 
tive relevance of these motions among the different samples 
of O VIII emitters. 

For this purpose we obtain the parameters /3 and g\2 
by minimizing the x^ASke function 
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Figure 9. Iso-corrclation contours for the measured two-point 
correlation function ^ in the plane (rp,7r) in redshift-space. The 
three panels refer to the three types of emitters considered. 
Bright-'WmM (top-left). Total gas {top-right). V\lYllM- dominated 
{bottom-left) The color scale quantifies the amplitude of the iso- 
correlation contours, {bottom-right). 



where 



y = log{l + £,{rp,Ti)), 



(9) 



Table 2. Best- fit values of the compression parameter /3 in the 
5(rp,7r) at rp > 6 Mpc h"! 



eij is the Poisson error in the measured correlation function, 
the upper script m indicates the model and the sum runs 
over all bins in which we have sampled the (rp, tt) plane. 

The distortion model provides a poor fit at small per- 
pendicular separation, as demonstrated by the fact that the 
reduced decreases significantly when one excludes pairs 
with rp < 6 Mpc h~^ from the analysis. This means that 
we cannot trust the value of (J12 obtained from the fit and 
therefore we cannot estimate the statistical significance of 
the Fingers of God-type distortions that appear more promi- 
nent in the WYiliA- dominated emitters. 

On the contrary, the model provides a good fit for 
rp > 6 Mpc h~^ and the best fitting value of /3 appears to 
be robust since it does not change when including/excluding 
pairs with rp < 6 Mpc h~^ from the fit or when we fix 
(T12 = 0. In Tab. [2] we list the best fit value of /? and its 1-cr 
error for the three types of emitters considered. We only con- 
sider the case in which we exclude pairs with rp < 6 Mpc h^^ 
since changing the cuts does not alter the result. For all 
types of emitters we measure a compression parameters P 
significantly above zero. This means that gas responsible for 
detectable O viii lines is typically found in regions that are 
moving coherently toward large density concentrations, as 
expected for O viii emitters that populate the outskirts of 
virialized structures rather than their central regions. Con- 
sidering the approximation involved in our error analysis 
and the uncertainties in the models of WHIM, we do not re- 



Emitter type 








Bright WHIM 


0.179 


± 


0.023 


Total gas 


0.152 


± 


0.018 


WHIM- dominated 


0.108 


± 


0.022 



gard the differences among the values of /3 for the different 
emitters as statistically significant. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have addressed two different, but tightly 
related, issues: the feasibility of indirect detection of the 
WHIM through the angular correlation analysis of the DXB 
and the study and characterization of the spatial distribu- 
tion of the WHIM already detected through its O viii X-ray 
line emission. 

In this paper we hav e use d the WHIM model described 
bv lBranchini et all (|2009l ') and lUrsino et all (|2010|) based on 
the hydrodynamical simulation of Borgani et al.l (|2004 ) with 
the metallicity-density relation taken from Cen & Ostriker 
(1999). Detectability and feasibility estimates assume next- 
generation X-ray satellite missions. As a reference case we 
have considered two instruments proposed in the framework 
of the proposed Xenia mission. For the angular correlation 
analysis of the DXB we have considered the HARI CCD- 
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detector, while for the spatial correlation study we have con- 
sidered the CRTS spectrometer. 

The main results of our analysis can be summarized as 
follows: 

• The DXB autocorrelation signal in the soft X-ray range 
[0.38-0.65] keV is mostly contributed by gas in hot and dense 
environments like clusters or groups. The WHIM contribu- 
tion is comparatively smaller but significantly different from 
zero. However, the intrinsic correlation signals of the WHIM 
and non- WHIM are similar, implying that the different con- 
tributions to the ACF of the DXB reflect the different SBs 
of the two components. No unique feature in the ACF of the 
WHIM can be used to unambiguously infer its presence from 
the measured ACF. The zero-crossing of the WHIM ACF 
at ~ 5' indicates the typical angular scale of the WHIM 
in emission and sets the minimum angular resolution re- 
quired to next generation instruments aimed at detecting 
the WHIM in emission. 

• According to our WHIM model, the non- WHIM sig- 
nal originates from compact regions associated to virialized 
structures, while the fainter WHIM emission traces the out- 
skirts of these regions. One can then hope to isolate the 
WHIM contribution by filtering out the signal associated 
to virialized structures, typically associated to the brightest 
spots in the maps. We found that the ACF measured after 
removing the 50% brightest pixels from the map is largely 
(~ 75%) contributed by the WHIM. More refined fihering 
procedures, in which one removes pixels associated to known 
objects, would certainly improve the results. 

• The presence of an isotropic and comparatively brighter 
Galactic foreground component represents the main chal- 
lenge to the WHIM indirect detection. We have found that 
a successful extraction of the WHIM ACF signal requires 
both an accurate modeling of the Galactic foreground and 
a long exposure time. The first requirement allows one to 
model the Galactic foreground as Poisson noise to be ac- 
counted for in the analysis. The second requirement allows 
to minimize shot noise errors. Our analysis suggests that 
exposure times of [O.f-f] Ms are required to measure the 
WHIM ACF. Uncertainties in the required exposure time 
generously allow for uncertainties in the WHIM model. Fur- 
thermore, it is possible to reduce the impact of the Galactic 
foreground studying the ACF at energies below ~ 570 eV 
(the local O vii triplet). This means focusing at emission 
lines at redshift significantly different from zero, thus re- 
moving the strong contribution of the Galactic O vii a nd 
O VIII foreground (as partly done in lGaleazzi et al^l2009^ . 

• The AGN contribution to the ACF is small and eas- 
ily removed by masking out pixels associated to resolved 
sources. In a deep (~ f Ms) exposure one expects to iden- 
tify ~ 2000 AGN in a f ° X f ° field. The corresponding frac- 
tion of removed pixels is a few % of the total. Unresolved 
AGN are basically uncorrelated and provides an isotropic 
contribution significantly smaller than that of the galactic 
foreground. 

• O VIII emitters that can be detected in a I Ms obser- 
vation show no density evolution in the range z £ [0,0.5]. 
Deviations from the constant density case reflect the increas- 
ing volume of the resolution element with the redshift and 
the patchy appearance of the "WRYiA- dominated component. 

• The two-point spatial correlation function of the O viii 



emitters can be measured with high statistical significance 
out to separations of ~ 10 Mpc h~^. Its amplitude and shape 
are close to those of optically selected galaxies, suggesting 
that both types of objects trace the same underlying mass 
distribution. O viii emitters contributed by the WHIM are 
more correlated, as expected in the framework of the peak- 
biasing model since they are preferentially found in the out- 
skirts of large virialized structures. 

• The correlation properties of O viii emitters do not sig- 
nificantly evolve in the interval z £ [0,0.5], despite the ap- 
preciable differential evolution in the mass fraction of the 
different gas phases. 

• Anisotropics in the clustering of O viii emitters induced 
by their peculiar velocities can be detected in the measured 
two-point correlation function and used to characterize their 
dynamical state. We have shown that the two-point correla- 
tion function of O viii emitters detected after an observation 
of f Ms will allow to unambiguously detect such distortions 
which, according to our predictions, should be characteristic 
of infall motions toward virialized regions. 

Our analysis indicates that next generation instruments 
will be able to detect the WHIM contribution to ACF of the 
DXB, but it will require a rather long exposure of ~ I Ms. 
This pr ediction seems t o be a t variance with the recent re- 
sults of iGalcazzi et"aLl (|2009l ') who found indirect evidence 
for a WHIM-like component in the angular correlation corre- 
lation properties of the DXB diffuse signal in several XMM- 
Newton fields with texp in the range [50-600] ks. As we ar- 
gued in Section [3.3.21 the discrepancy is probably apparent 
since our predictions are to be regarded as conservative for 
two reasons. First of all, our model overestimates the SB 
of the non- WHIM components hence artificially increasing 
their contribution to the total ACF compared to that of the 
WHIM. Second of all, we use a very simplistic cleaning pro- 
cedure that is not optimized to extract the WHIM signal. A 
better treatment of the non- WHIM component, capable of 
bringing theoretical prediction into agreement with currents 
DXB constraints, and a more effective cleaning procedure, 
that pre ferentially removes p ixels associated to known ob- 
jects, as iGaleazzi et al.l (|2009l ') did, will significantly reduce 
exposure time required to per form the angular cor relation 
analysis. With this respect the IGaleazzi et al] (|2009l ) results 
give us confidence that next generation instruments will be 
able to carry out a very effective observational campaign for 
the indirect search of the WHIM within a reasonably limited 
amount of time. 

In fact, the comparison with IGaleazzi et al.l (|2009l ) re- 
sults constitutes a new obser vational test fo r our WHIM 
model. The ACF measured bv IGaleazzi et all (|2009l 'l (their 
Fig. f ) is in agreement with the one predicted by our model 
(Fig. O for angular separations larger than 9 = 2'. At 
smaller separation the mismatch simply refiects the smooth- 
ing effect of the instrumental PSF which we did not try to 
model in our mock maps. 

The fact that our WHIM model passes this further ob- 
servational tests(the others being the O vi and O vii ab- 
sorbers abundance and the DXB intensity) is certainly reas- 
suring, but is no proof that our model is correct and its pre- 
dictions robust. This concern, common to all WHIM mod- 
els, has triggered a number of works aimed at testing the 
sensitivity of model predictions on crucial, ill-known pro- 
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cesse s of metal diffusion and stellar feedback (iBertone et al.l 
I2OI0I : iTornatore et aljboiol : IShen et al.ll2010l ). Focusing on 
emission, we have seen that the main inconsistency of our 
model is the excess SB of the non- WHIM gas, whereas the 
soft X-ray emission of the WHIM meets the observational 
constraints. The natural question, then, is whether changing 
theoretical prescriptions for metal diffusion and galaxy feed- 
back may help reducing the non- WHIM emission. One pos- 
sibility is to invoke an AGN- driven rather tha n supe rnovae- 
driven feedback mechanism. ITornatore et al.l (|2010l ) found 
that, on average, AGN feedback increases the mass frac- 
tion and metallicity of the WHIM in regions with overden- 
sities between a few and 10. However, the bulk of detectable 
WHIM emission is produced in regions with higher overden- 
sity, in which AGN feedback decreases both the mass frac- 
tion a nd metallicity of the gas. The analysis of lBertone et al.l 
is perhaps more relevant to us, since it focuses on 
the O VIII emission line. They show that detectability pre- 
dictions are quite robust to the different model prescrip- 
tion with the exception of the cooling function. They show 
that using a metal-dependent cooling instead of a metal- 
independent scheme (as in our model) decreases by up to an 
order of magnitude the SB of the O viii line in the bright 
areas {SB ^ 1 photon cm~^ s~^ sr"'^), dominated by non- 
WHIM emission, but does not significantly alter the SB of 
the fainter areas, largely contributed by WHIM emission. 
Then, on a qualitative level, the use of a metal-dependent 
cooling scheme, would help relieving the problem of our 
models. 

Assessing, in a quantitative way, the impact of these re- 
sults on the possibility of detecting and studying the WHIM 
is beyond the scope of this paper. Here we have simply es- 
timated their effect by modifying our WHIM model in a 
way that effectively (but qualitatively!) mimics the impact 
of using a metal-dependent cooling scheme. For this purpose 
we have assumed an ad-hoc density-metallicity relation such 
that the metallicity of the gas with WHIM density remains 
unchanged, while it is decreased in high density region. We 
stress the fact that this is an ad-hoc model that only serves 
for illustrative purposes. The result of using this metallicity 
scheme is to decrease the mean SB in high density regions 
by ~ 30% while the one of the WHIM remains the same. 
The net effect is to reduce the SB contrast of the two phases 
although the relative statistics does not change appreciably, 
as indicated by the fact that in both the new and origi- 
nal metallicity schemes 50% of the pixels are dominated by 
WHIM emission. In terms of correlation analysis, the angu- 
lar ACF of the WHIM remains the same but its statistical 
significance after the cleaning procedure increases signifi- 
cantly, as expected. 

The results of this work and the robustness of the 
WHIM model predictions, suggest that a next-generation 
mission like Xema will be capable of extracting and charac- 
terizing the WHIM contribution to the angular correlation 
properties DXB. We also expect to detect enough O viii (and 
O vii) emitters to characterize the spatial correlation prop- 
erties of the WHIM and its dynamical state. We confirm 
that, even with long exposure, next-generation instruments 
will only be capable of studying the spatial distribution of 
the WHIM in metal-rich, high density regions associated to 
the outskirts of large virialized structures, and of detecting 
the signature of its infall motion. Baryons tracing the large 



scale filaments of the cosmic web will remain largely unde- 
tected. 
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